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Abstract 

A new scheme has been proposed to solve the B. E. condenstates in 
terms of Green's function approach. It has been shown that the radial 
wave function of two interacting atoms, moving in a common har- 
monic oscillator potential modified by an effctive interaction, satisfies 
an intregral equation whose kernel is separable. The solution of the 
integral equation can be written in terms of the harmonic oscillator 
wave functions. The ground state wave function of the system can be 
written in terms of these solutions. 

PACS numbers: 03.75Fi,02.60Cb,32.80Pj Keywords : B. E. Condensate, Ef- 
fective Interaction, Green 'Function. 
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The recent experimental realization of Bose-Einstein condensation [BEC] 
in ultra cold atomic gases has generated intensive theoretical and practical in- 
terest regarding the properties of Bose gas. Non-linear Schrodinger equation 
in this context known as Gross-Pitaevskii [ G-P ] equation [l],has been used 
to study the properties of the interacting spinless Bosons , trapped in har- 
monic potential . Altough the G-P equation is widely accepeted as a workable 
model for description of the B E condenstates at low temperature, the actual 
dynamics of B E condenstates is not fully captured by the G-P equation be- 
cause of the non-integrability of the G-P equation and explicit solution is not 
known [2]. It is to be noted that for the positive scattering length, the G-P 
equation is generally solved numerically only under cylindrically symmetric 
systems. The Thomas Fermi approximation [3] has been introduced to get 
some analytical solution for the description of B E condenstates. On the 
other hand for the negative scattering length the G-P equation is yet to be 
properly investigated [3]properly. Another alternative approach to study the 
dynamics of the B.E. condenstates is the time dependent variation technique 
[4] in which one assumes fixed profile of the condenstate and computes the 
evolution of some parameters such as the width by variational techniques. It 
should also be mentioned that the negative scattering length provides stable 
solution of the G-P equation only under certain conditions for the number of 
particles and the size of the trap [5]. When those conditions are not fulfiled 
the condenstates become unstable. Therefore realization of large conden- 
states will provide reliable experimental data for the precise understanding 
of the condensation process[6]. In this paper we shall present a new solution 
for the condenstate on the basis of Green's function approach of scattering 
theory. Similar method was used in the treatment of the nuclear many body 
problems by Bethe and Goldstone; Brueckner and Gammel [7] for the un- 
derstanding of the various aspects nuclear many body problem. We have 
derived an integral equation for the description of two interacting bosons in 
the presence of many bosons. We shall consider a system of N spinless in- 
teracting atoms all having the same hyperfine species, in equilibrium and at 
very low temperature and densities characteristic of alkali atoms under Bose- 
Einstein condensation (BEC) condition. We assume that each of the atoms 
is moving in the influence of a common harmonic oscillator potential which is 
called a trapping potentia^V^) and an effective two-body interaction acting 
between them which modifies their behaviour . The essential questions we 
wish to investigate are : How does the effective two body interaction modify 
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the two-body wave function?. How does the energies of the two interacting 
bosons system change when the two- body effective interaction is switched 
on ? It is welknown that the scattering length (a s ) is very small compared 
to the interatomic distant (r int ) [8] under BEC in the dilute atomic alkali 
gases. Therefore the behaviour of the atoms is extremely sensitive to the 
details of the trap and the effective two- body interaction. Since the effec- 
tive two- body interaction is very weak the quantitative calculation based on 
perturbation theory in the interatomic interaction is hoped to be highly reli- 
able for the alkali gases. Furthermore the de Broglie wave length associated 
with individual atom is comparable to interatomic Ti nt distance therefore the 
correlation between atoms becomes important [7]. Hence the construction of 
the many -body condenstate state wave function of the system in terms of 
two-body wave function will be interesting. The standard many-body wave 
function in terms of Hartree-Fock ansatz using the single particle solution of 
the G-P equation does not take into accoutn the effect of correlation between 
atoms which is very important in the case of BEC alkali gases. In this letter 
we have presented an analytical formalism to investigate the above aspects 
of BEC problem without any numerical estimate for a particular choice of 
the effective interaction. That will be communicated shortly. 

Let us consider a system of interacting N spinless atoms and now arbi- 
trarily pick out any two interacting spinless atoms, say 1 and 2, denoting 
their relative and center-of- mass co-ordinates, by f and R respectively . We 
shall investigate the state (0) of these two interacting atoms trapped in a 
common harmonic potential in the presence (N — 2) interacting atoms. The 
most general states of the many-body system will then be described by a 
wave function \l/(r, R, x, t) where x represents schemitically the co-ordinates 
of the remaining (N — 2) particles. It is the dependence on r that is crucial 
but the parameter dependence of \& on R, x and t will be ignored . In the 
present context we assume that simple BEC is realised in the two-particle 
state (f>{r) whose equatiom of motion will be determined using Schrodinger 
equation. From the knowledge of the wave-function of two-atom system (f>(r), 
the many- body ground state of B E condenstate system of alkali atoms can 
be written. The hamiltonian for the two- atom system in the common har- 
monic oscillator potential of frequency co is given by 
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Introducing the cannonical co-ordinate transform 



and 



r = r 2 -ri,p= , (2) 



R= { ^p^-,P = (P2 + Pi). (3) 

The hamiltonian is transformed to 

(P 2 + M 2 ^ 2 ) ( P 2 + M W) 
^ = 2M + 2^ ' (4) 

where M=2m,^ = y. The Schrodinger equation for the system with % — 1 
is 

[ ( -v| + mw> + ( -v; +^v) 10(r>( g) = ^3 (5) 

The Schrodinger equation for the relative motion is then 

= {~ + kr 2 )4>(r) = E ^(r), (6) 

where 

* = ^ = ^. (7) 
The solutions are given by 

Rniir 



if) = -^^Y lm (r), (8) 



where Yi m (r) are the usual spherical harmonics. The radial wave function 
R n i(r) is given by 
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Rm{r) = N nl exp{- — )r l+1 v nl {r), (9) 
where v — ^ and f„;(r) is the associated Laguerre polynomial: 

Mo=4 + 4 +i K 2 ). (io) 



N n i is determind by the normalization condition and E^i = hu(2n + 1 + |). 

We shall now treat the effective spin independent two-body interaction 
v(r) which is very weak between two atoms , as a peturbation and solve the 
Schrodinger equation using the Green's function method to determine the 
perturbed wave function for the two- atom system. The unperturbed wave 
function is given by 

Mrl = ^-Y lm (r). (11) 
r 

Similarly the perturbed wave function can be written as 

^ nih (r) = ^Y hmi (f). (12) 
r 

The Schrodinger equation for the perturbed system is 

[H Q + v{r)p±Y mih {r) = E nih ^Y mih (f) , (13) 
where E mil is the perturbed energy eigen value. We now write 

H Q ^Y lm (r) = ~[-~ + V^RUrWUr) = E° nl ^Y nl (r). (14) 
r r m or r 

Multiplying the perturbed equation(13) by Y* ^ (f) on the left and integrating 
over the solid angle dr we obtain 

1 d 2 

[Ernh ~ + ^i( r )Ki«i( r ) = v h (r)u nih (r), (15) 

where 

v h (r) = J drY* mi (r)v(r)Y hmi (r) (16) 

In order to solve this equation we introduce Green's G ni i'(f, r 1 ) corresponding 
to equ.(15) satisfies the equation 

(E nih - H )G nih (f,f) = 5{r - f). (17) 

Multiplying through on the left by ^* mi (r) and integrate over the solid angle 
dr we obtain 

[E nih ~ (-^) + V h (r)}g nill (r,r') = S(r - r'), (18) 



where 

g nil Ary) = Y, Rn2l i r)Rn2h{r '\ (19) 

where 

E mli,nih = (E ni h ~ E n 2 h)- ( 20 ) 

Hence the perturbed wave function u niJl (r) satisfies an integral equation 
with a separable kernel. Therefore we write 

roo 

M m/i( r ) = / dr'g nih (r,r')v h (r')u nih (r'). (21) 
•/ o 

Finally we write 

UnM = £ * n2h(r ) rdr'R n2i y)v h {r')u nii y). (22) 
n 2 =o E ni h,n 2 h Jo 

The equation for R nll (r) is 

[^ 1 -(- 1 ^ + ^ 1 (r))]i? nil (r) = o. (23) 

Multiplying equ.(15)and equ.(23) on the left by i? n ; x (r) and w niil respectively 
subtracting and integrating over r' we obtain 



as 



and the perturbed wavefunction « ni ^(r) in a closed form can be written 

oo 

u nih( r )= Yl K nih, n 2 liRn2h ( r ), ( 25 ) 
n 2 =0 

where 

1 r°° 

K ni i un2h = — / drR n2h (r)v h (r)u nih (r). (26) 

*- J nih,n2,h •'O 

. To this expression of we should add the solution of (-Ey^ —Ho)Xn 1 i 1 = 
0. Since i?™^ 7^ the equation (E ni ^ —H )x ni h = has no solution. Thus 
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it is not necessary to add a particular solution to the . From equa- 
tions^) and (26)it follows that 

K ni h,n 2 h = p J2 K nill , nsh / drR n2h (r)v h {r)R n . ih (r). (27) 

- C/ raiii,n 2 /i n 3 =0 

Results and Discussion 

The main conclusions of the paper are the following: It should be noted 
that 77-3 can in principle run from o to oo, the above equ.(27) reprsents an 
infinite set of simultaneous homogeneous algebraic equations in K n] i x<n ^. 
The consistency of the these equations imposes the condition of vanishing 
of the determinant formed from the coefficients of iT ni ; i n2 ; 1 . This condition 
determines the energy shift E n ^ n \^oi the two interacting atoms. The corre- 
sponding eigenvectors when substituted in eq.(26) determine the perturbed 
radial function u ni ^(r). Since B E condensation takes place strictly in the 
ground state of the system therefore a few values of n 3 in the sum will pro- 
vide a very reliable estimate of the two-particle wave function. Because of 
the separable structure of eq.(21) we can also include various types of central 
and non central of two-body effective interaction and the variation of the 
scattering length can also be incorporated to investigate sensitivity of the 
two-body condenstate wave function and the stability of the Bose-Einstein 
condensation in alkali atoms on the effective two-body interaction which can 
not be so easily taken care in the formalism of G-P equation. When the 
effective two-body interaction becomes non-central in nature then the radial 
equation for u ni i x {r) i.e eq.(21) will satisfy a coupled integral equation. When 
the effective two-body interaction is a pseudo- potential type, as commonly 
used in the calculation of G-P equation, the determination of the two-body 
wave function Un^ir) becomes extermely simple in this formalism. We can 
also investigate the importance two- particle correlation in BEG Finally the 
construction of the condenstate wave function describing the BEC of the in 
terms of the calculated two- particle wave function to study various aspects 
of BEC becomes possible. These are the main conclusions of the paper. 
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